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MAGNITUDE OF NONLINEARITIES IN COLLISIONLESS MAGNETOPLASMA WAVES 

by Richard R. Woollett 
Lewis Research Center 

SUMMARY 

A perturbation technique is used to obtain approximate solutions to the nonlinear 
hydromagnetic equations (Ohm's law and the equation of motion of a fluid element) of a 
cold collisionless plasma. The assumed model consequently consists of (1) a linear 
harmonic wave that is infinite in extent and (2) nonlinear waves that are harmonic and 
propagate in the direction of the fundamental. The perturbation technique limits appli- 
cation to regions where nonlinearity is small. Some idea of the importance of the non- 
linearity can be obtained by observing the magnitude of the terms. For a transverse 

perpendicularly propagating wave in atomic hydrogen, nonlinear effects are shown to be 

7 -1 

important near and below a frequency of 10 second for a fundamental amplitude of 

12 

1 statvolt per centimeter or greater and for a plasma density of 10 ions per cubic 

4 

centimeter and a magnetic field intensity of 10 gauss. For a longitudinal perpendicularly 
propagating or a right or left circularly polarized wave, nonlinear effects are only 
important for very large wave amplitudes. 


INTRODUCTION 

Nonlinear calculations for plasma wave propagation are generally very difficult be- 
cause of the complexity of the hydromagnetic equation (Ohm's law and the equation of 
motion of a fluid element). Consequently, to make the problem tractable, simplifying 
assumptions must be made. One such heuristic hypothesis entails assuming a specific 
structure for the nonlinear plasma wave. Experimental characteristics can be used as 
a guide to establish the a priori structure. One typical characteristic of experimental 
nonlinear steady -state systems is that numerous discrete frequencies in addition to the 
frequency of the forcing oscillation (pp. 81 and 82, ref. 1) are frequently present. With 
this information as a guide, the following model was constructed and used herein. 

The linear wave present is assumed to be a single plane wave. The nonlinear dis- 
turbance is assumed to consist of an infinite set of plane waves, which are periodic in 


both space and time and propagate in the direction of the linear component. The order 
of magnitude of the ascending coefficients of the series is monotonically decreasing. A 
perturbation technique (ref. 1) that can readily incorporate these assumptions into the 
nonlinear equations is used. With the model so defined, it is possible to calculate the 
relative magnitude of the nonlinear effects. 

The theory is developed for plasma waves in an unbounded medium propagating both 
perpendicular and parallel to a superimposed steady -state magnetic field. 

A similar analysis has been made (ref. 2) for waves propagating perpendicular to a 
magnetic field immersed in an electron gas that has a uniform neutralizing background of 
infinitely heavy ions. In the present calculation, however, the effects of the cooperative 
motion of singularly charged ions are included; in addition, the case of waves propaga- 
ting parallel along the magnetic field is investigated. Calculations of the amplitudes of 

the first perturbation of the linearized wave in atomic hydrogen are presented for 

2 3 4 

frequencies near the cyclotron conditions for magnetic fields of 10 , 10 , and 10 gauss 

11 12 14 _3 

and plasma densities of 10 , 10 , and 10 centimeter . The results are also pre- 

sented in nondimensional form. 
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Subscripts: 


electron 
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x,y,z 

Cartesian components 

0 

steady -state conditions 

1,2, 

order of perturbation 

and 3 
1 

perpendicular to magnetic field 

II 

parallel to magnetic field 


Superscripts: 

- vector quantity 

~ factor of quantity independent of space and time 

~ unit vector 


THEORY 


The major characteristics of plasma waves can be determined from the dispersion 
relation that relates the wave number to the forcing frequency. The dispersion relation 
of a magnetoplasma can be derived with the use of Ohm's law and the equation of motion 
for a fluid element of ionized gas. These two expressions, however, are generally not 
referenced in a convenient form that includes the nonlinear terms. In order to obtain 
these terms, Ohm's law and the equation of motion are derived starting with the hydro- 
magnetic equations (ref. 3, eq. 2-4) of the various charged species present in the 
plasma. For the ions, this expression is 


l dv i _ _ \ h v i xB i — - 

n l m i( — + v i ' Vv i) - n i e ( E + -7- ) - v ' *1 * p e, i 


( 1 ) 


and for the electrons, it is 


'3v 


v„ x B 


VM- if ‘ = ' n eT + c 


- V '*e + P e,i 


( 2 ) 


Adding equations (1) and (2) while making use of the conservation of mass 
dp/dt + V ’ pv = 0 results in the following form of the equation of motion: 


6 



( 3 ) 


p $1 + pv ’ Vv + (if — (J • VJ + JV * J) = J X B + P e E 
3t \e / n 


When equations (1) and (2) are multiplied by e/cm. and e/cm. and are then subtracted, 
the following form of Ohm’s law results: 


me -7 m a 
e 3J . e 


m_\ in 
e\» . e e — 


e 


2 at 


+ — v ’ < n i v i v i - n e v e v e> = n e 1 + ~ E + E + 


vx B - ill 


m. 


m. e 


m^c 



x B 


where v and J are defined by 

pv = n e m e v e + n-m-v. 


(4) 

(5) 


and 


; , = n iVV e (6) 

e 

Since the plasma is a cold and collisionless gas, the stress tensor and momentum trans- 
fer term have been neglected in the derivation of equations (3) and (4). It is convenient to 
derive the nonlinear effects for the perpendicular propagating wave (k perpendicular to 
Bq) separately from those of the parallel propagating wave (k parallel to Bq). 


Perpendicular Propagation (k Perpendicular to Bq) 

To specify the form of the solution, it is first assumed that the running wave of the 
fundamental along with any higher or lower harmonics all propagate perpendicular to the 
magnetic field; that is, k perpendicular to Bq. Further, it is assumed that any forced 
disturbance shall be composed of traveling waves of the form 


A = 




i(k.x+w-t) 

J J 


(7) 


where j = 1 is the linear solution of equations (3) and (4). The components of the series 


7 



for j >1 originate from the nonlinear terms in equations (3) and (4). If an additional 
assumption is made stating that terms in ascending order of j represent terms of de- 
creasing orders of magnitude, it becomes possible to form a set of recurrence equations. 
If the zero-order values of p and B are constant while those of J, p g , v, and E 
vanish, the first-order terms of equations (3) and (4) reduce to 

3vi - - 

Pn = Ji X Bq (8) 

u at 1 u 


and 


m e v 1 x B 0 , _ _ 

— = E-, +— - — J, x B n 

2 at 1 c ne 1 U 
e n 


(9) 


respectively. 

The assumption of an exponential solution (eq. (7)) permits gradients and time 
derivatives to be expressed as ikx and iu>, respectively, so that v • A and 3A/at 
now becomes 


OO 


OO 


2>i 


A. and 
h A 


EVr 


respectively. 


These equivalences are used throughout the 


j=l 3=1 

the following derivations. From the conservation of mass, 


Pi = -P 0 



l,x 


(10) 


Consequently, the second-order component of equation (3) is 


3v 9 k« 

Pn — - Pn — v 


9v, 


at 


+ Pn v i • Vv, + 


vvi 


/ \2 m„ 


VJj + JjV • Jj) 


- ^2 * ^1 * Pi, e ®i (H) 


The charge density is expressed as 


p l,e = e < n l,i - n l,e> 


( 12 ) 


8 



but since 


ki m„ „ 

ri — 1 „ T - G C » 

n l i — — n 0 V 1 x + "lx 

1,1 u i,x p e ’ 


(13) 


and 


k l m i c 

n l e = " — n 0 V 1 x “ J 1 x 

A,e u 1,x p e ’ 


(14) 


which are obtained from the conservation of the mass of the species, p. may be 

1 , G 

written as 


Oi,. 


(15) 


Substituting equation (15) along with the identities 


V 1 • Vv l = v l,x ik l v l 


(16) 


and 


V vjj ♦ J,V • Jj - 


(17) 


into equation (11) results in 


c 2 m k, _ i _ _ k i p - 

v« = ~2 — J J - _i_ J,xB n - -1 l — x B« - — t £- J, E 1 

e n Q p 0 a . 2 X ’ x P 0 u > 2 P 0 " 2 P 0 W 2 iw l 


(18) 


When equations (14) and (15) are used, the second-order equation obtained from Ohm's 
law reduces to 


9 



m 


-1T + T- V ' - Wl,e> = E 2 - - (v M 

n a i 


m e k l c . - v, x Bj v 2 x E 0 

J i Cj -j + — 4- 

m i w i n o e ’ c c 


k„v 


'4rr v i XB o-ri J i XB i- J - j 2 xb o < 19 > 




n O e 


n O e 


The spatial derivative term of equation (19) can be replaced by 


v ’ (v l,i v l,i- v l,e v l,e ) - 


2ick 


n Q e 


“( v l J l,x + J l v l,x 


m.c 

Pe 


J l,X J l y 


( 20 ) 


Expressions for^Vj and Vg are given by equations (8) and (18), respectively. From 
Faraday's law, Bj can be represented as 


_ ck., 

B. = — t (E, y - E. z) 
1 w v 1 , l,y ; 


( 21 ) 


When these substitutions are made, equation (19) becomes 


iB o _ 


lV2 'iw J2,i ' E2 


m 210^/180 . _ »B 0 m.c _\ ik.B 0 _ k.m.c 

= — — (Jj X y + J t yX)J x +— “ Jj J ) +— i-5 J. E, + JLi.j 

e n Q e \p 0 «l ’ ’ r ’ <>0 w l ’ V * *’* 7 PfflW J hV 1 ^ h 


x E l 


2cm Q k 1 B fl r A 

W^l.lEl.y^l.yEl.z 2 ) ' [hj'hy + J l,y x >] 


Vl 


"O^O 


lk l B 0 E 1 , y j . lk . l E l,z J l,z B 0 > + lk l J l,x B 0 ( E J + E j) . k . l E 0 J l»y J 

’ Po u 2 w l v 0 u>iP n 1,x 1,y , ,3„„2 1 ’ 1 


2“ 1^0 


w l cp 0 


1 r ~i Bn 

— L E l,y( J l )X ^ -Jl.^ + El.z^l.x 2 - J l,z«] - — (- J 2,xy + J 2,y*) (22) 


“0 C X 
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In order to compare the amplitudes of the nonlinear fields with those of the linear fields, 
it is necessary to first express Jg as a function of and E 2 * This operation is 
accomplished by substituting the following expressions for Jj (obtained from eqs. (8) 
and (9)) into equation (22): 


io>. 


E l,x n i " E l, y 


a co 


l,x 


C,1 


E 1 + E - 

l,y I l,x 


xco 1 


au ) 


i>y 


C,1 


j 


1, z 


u> 


cia>. 


17 


(23a) 


(23b) 


(23c) 


Before proceeding, it will be convenient at this point in the derivation to split the 
problem into two cases. These two cases arise since there are two dispersion relations 
for the linear wave equation. The wave equation can be expressed as 



Since is a function of E<, equation (24) represents a set of three homogeneous 
equations. For perpendicular propagating waves (k perpendicular to Bq), Ej is not 
coupled to E t or E, ... Consequently, two independent forms of the dispersion 
relation exist. The wave consisting only of E, is called the ordinary wave, while the 
remaining wave consisting only of E ± is the extraordinary wave. The dispersion rela- 


11 


tion for the ordinary wave is 



( 25 ) 


while for the extraordinary wave it is 


k 


2 

1 



(26) 


These expressions are identical to those presented in reference 4. Since both waves are 
independent, the higher-order effects will be developed first for the extraordinary wave 
and then for the ordinary wave. 

Extraordinary wave (k perpendicular to B Q , Ej perpendicular to B Q ). - If the 

linear currents expressed by equation (23) are substituted into equation (22), the follow- 
ing expressions can be obtained for the first harmonic currents when Ej z = 0: 


fi II E 2,x 


lQJr 


aw 


E 2,y + S(J l,x’ J l,y ) 


%x 


c,i 


(27a) 


lW f 


'2,y 


% E 2,y + E 2,x + T(J l,x’ J l,y ) 

c,i 




2, z “ t : ^2, z 

’ l47TCO)o 9 


(27b) 


(27c) 


where 


® s x J l,x + *x**l,y + “x^l.y^ljX + x x^l,x + ^x^l,y 


(28) 


12 



and 


(29) 


T " s y J l,x + V J l»y + u y J l,y J l,x + x y J l,x + y y J l,y 

- i^x+Wot) 

Substituting the expressions for J 2 into the wave equation of the A 2 e wave 

gives 


a) 



2,x 


-1477 


Cl), 


'2,y 


= 0 (30) 


*2,y 


This substitution results in a set of nonhomogene ous equations that can be solved for the 
electric fields. Hence, 


E 2,x “ E 2,x e 


i^x+ajgt) 






i(2kjX+2u)jt) 

e 



(31a) 
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E 2,y = E 2,y e 


i(k 2 x+w 2 t) = m (<*“ c,i 


IV 


u>. 


fi m E 2,x + ~ ) e 


g \ i(2k 1 x+2u> 1 t) 


( 3 * 


IV/ 


f ^- k 2 !ilE n _ = 0 


2, z 


(3 


2 -1 
Both S/Ej x and T/Ej x are proportional to Bq ; in addition, they are functions of 


Wp Wp, and W g only, where 


W = 


Cl) 


U) 


c, i 


Cl) 

w =-® i £ 

e U) . 
C,1 


0) 

W =_P 
p 


Cl) 

C,1 


Consequently, these terms can be written as 

— S(W,W ,WJ 
£2 B n P e 


and 



= — T(W,W ,W ) 
Bn P 6 


Equati° ns (31a) and (31b) can thus be nondimensionalized when multiplying through by 
B n /E ? , which results in 



and 


E 2,x B 0 _ 


SI 


IV 


S(fi m - Kj) + IT 


W 2 W e 

w p n ivj 


E 


l,x 


? / w P W p 

- K V - 1 , 

w pV 


(32a) 


E 2, y B 0 _ _ ^IV W p 
g? " ' W e W 2 

1 , X 



2,x B 0 


1, 



(32b) 


where K(= kc/w), fi m , and fl TV are nondimensional expressions and functions of W, 

ill ro 1 V ^vn /v (vO 

W p , and W e only. Since E 2 x Bq/E^ x and E 2 y B 0 /Ep x are invar ient to any changes 
in Wp n, Bq, and m i where Wp W p , and W g remain constant, they are convenient 
forms in which to present the results. 

These expressions for the various second-order electric fields are functions of x 
and t. Since the amplitudes of the various harmonics are constant, the only way that 

t is for 


(33) 


E 2 and Ep Substituting these 

E 0 must vanish. Another 
Z j z 

represented by E 0 „ and E 0 „ 

E 9 can be determined from 
(31c). However, the resulting dispersion relation is just that for the ordinary wave ob- 
tained from linearized theory. Since it is not coupled to the nonlinear wave and since it 
is specified that the fundamental is an extraordinary wave, E 2 z can again be set equal 
to zero. The electric fields associated with the second-order perturbation, consequently, 
are perpendicular to Bq and are functions of the superimposed steady-state magnetic 
field. In addition, since all disturbances propagate transverse to the magnetic field, 
the first harmonic may be called an extraordinary wave. 

In summary, when an extraordinary wave exists in a plasma, the nonlinearity of the 
equations requires an additional extraordinary wave. This nonlinear component has the 


equations (31a) and (31b) may be satisfied for all x ar 


k 2 = 2kj 


u> 2 = 2o>i 


The phase velocity, consequently, is identical for both 

values of 1^ and a> 2 into equation (31c) dictates that 

possible solution exists for equation (31c) . If the wave 

is uncoupled from E 0 _, a separate wave number for 

z * z 
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same phase velocity but twice the frequency as the linear component. 
The wave amplitude, that is, 


>eE 2x ) 2 + «ReE 2v ) 2 ] 


1/2 


2,y' 


or 



(34) 


is probably of more interest than the expressions for the nonlinear components of the 

disturbance. This equation is still inconvenient because of the terms E, and E, 

a > x a f y 

However, E, ,, is directly proportional to E, v (see eqs. (23) and (24)) so that a 
9 a > y A, X 

factor of Ef may be extracted from equations (31a) and (31b) and hence cancelled 

A j X 

from the resulting expression. If the right hand side of equation (34) is divided by 


E 


l,x 



which is proportional to 



2 


there results 



(35) 


(36) 


which is an expression independent of Ej x and E^ y. Both j El 2 1 and |Ej| repre- 
sent the time average of the wave amplitude. Equation (36) will hereinafter be referred to 
as the ’’amplitude ratio” of the extraordinary wave. It is the ratio of the amplitude of 
the nonlinear wave to the square of the amplitude of the linear wave. A nondimensional 
wave amplitude ratio | E 2 1 Bq/ |Ej|^ will also be investigated. The description of the 
behavior of the amplitude ratios will be discussed in the section RESULTS. 

Ordinary wave (k perpendicular to Bq, Ej parallel to Bq). - In this case, 

Ej x = Ej y = 0 while Ej z * 0. In an approach similar to that used with the extra- 
ordinary wave, plasma currents may be determined for the ordinary waves. Thus, for 

E l,y = E l,x = 0 and E l,z^ 0 - 


16 
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(38c) 




= 0 


~o 

Equations (38a) and (38b) can be nondimensionalized by multiplying through by Bq/Ej z 
and by expressing the variables in terms of Wp Wp, and W g . This results in 


E 2, x B 0 _ 


K l W e 

n iV W l W p 


W^SJjji - K 2 2 ) , 


w. 




(W| - 1) 


(39a) 


E 


l,z 




m 


(%I - K 2) 


¥= 

‘^IV, 


and 


E 2, y B 0 _ 


E 


l,z 


K l W e 

n lV W l W ? 


a 


m (w 2 , 


w 0 


D + 


w; 


w: 


pJ 


(39b) 


l 2\ l W 2 W e 

%i( n m - K l) J 


w p"iv/ 


(VI (V <\ I*V (V/O 

Since E 2 x Bq/Ej z and E 2 yBQ/Ej z are invarient to any change in Wj, n, B Q , and 
m i where Wp W p , and W g remain constant, they represent a convenient form in 
which to present the results. 

These second-order electric fields are analogous to the fields of the extraordinary 
fundamental, since both depend on the magnetic field and since both have x and y 
electric field components but no z component. Consequently, the first harmonic 
associated with an ordinary wave fundamental can be called an extraordinary wave. The 
propagation constant k 2 and the frequency w 2 can be shown to be equal to 2kj and 
2u>i in a fashion similar to that used with the fundamental extraordinary plasma wave. 
Substituting these values of 1^ and o> 2 into equation (38c) dictates that a coupled E 2 z 
wave must vanish. An uncoupled wave E 0 cannot exist for the same reason as was 
presented in the previous section. The amplitude ratio for the ordinary wave, conse- 
quently, is 


18 



( 40 ) 



I /~ i . />-> | o 

E 9 / E, „ . A nondimensional wave amplitude ratio 

/v ~ n « 1 ) z 

I E 2 1 Bq/ I E j p will also be investigated. The description of the behavior of the amplitude 
ratios will be discussed in the section RESULTS. 


Parallel Propagation (k Parallel to Bq, E 1 Perpendicular to Bq) 

When k is parallel to Bq, equation (7) is not valid; in this case, 



~ i( k i „z+w,t) 


V 


]>z 1 


(41) 


Since the propagation vector k is now directed axially, all relations involving k and 
its spatial derivative must be changed. The derivation of the amplitude ratio is presented 
with equations analogous to those used with perpendicularly propagating waves. With the 
use of 


v ' v l= ik l,z v l,z 


the zero-order equation of motion becomes 


P 1 = " p 0 


Ail 

Wi 


1,Z 


(42) 


while the second -order part of the equation of motion becomes 


P 0 1W 2 V 2 + Dlk l,z J l,z J l - J 2 x B 0 + J 1 x B 1 ‘ ~~ cJ l,z E l 

W 1 


(43) 


When equation (12) and 
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k l z n 0 / m e c 

"i,r-^Kz*7 J U 


k l,z n O L m ic. 

V»-*'7e M 


are used, the charge density may be written as 


p .’ e = '7 icJ i’ z 


Substituting equation (46) along with 


V 1 ■ Vv l' lk l,z v l,z v l 


J 1 • vjj - Jj v • - lk lf z J 1>z Jj 


into equation (43) results in 


- _ " Dk l,z J l,z J l iJ 2 XB 0 iJ i x B i ik i,z cJ l,z E l 
v o “ 


P 0 W 2 


p 0 W 2 


P 0 W 2 


p 0 w l w 2 


When equations (45) and (46) are used, the second-order equation obtained from Ohm's 
law reduces to 


m e 0J« m / v _ L , / m. \_ 

+ ~ v ‘ ( v i i v i i " v i o v i = E 9 - (v- - __L cJ, )e , 

J n 9t e V 1,1 lfl 1>e 1 » e ' 2 cdj \ 1,z p Q e 1,z / 1 


m e k l,z c T - . v l XB l . v 2 xB 0 p O k l,z v 


m i n 0 ew l 


v r0j v 2 , D() h 0 1, z 1, z — 
l,z E l + + ‘ V 1 


n O m i cco l 


n 0 e n 0 e 
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The spatial derivative term can be replaced by 


2icm.k« 


- m.c _ 
l,z v l + v l,z J l ■ ~ J l,Z J lj 


( 51 ) 


Expressions for Vj and v 2 are given in equations (8) and (49), respectively. From 
Faraday's law, Bj can be represented as 

— ck i „ 

B l=-7 £ (El >y *- E l > ^ ) < 52 > 

W 1 

When these substitutions are made, equation (50) becomes 


im^cw, _ m /2ick, \( m.c 

— - — — Jo + — -I ^IIJ, ,Vi „j,l= E 


e^n. 


n o e 


X ’ Z 1,2 ir 2 


ck, m, _ v.xB, Dk, B n 

+ h z Ej ♦ -i— i - hi 0 j (_J y + J v) 

i-o/.i A > z 1 A ’ z A ’ x A >y 


p o ea, l 


cp n w 


0 2 


iBj _ ik, _B, 


cp 0 w 2 


— AXV 1 A A 

J 2 , - -- - (-Ji z Ej x y + J x Z E X x) 

2,1 pQWgWj 1,2 1,x 1,2 1,r 


ik-t „ B A A A ^ J 2 ^ Bq 


axv i 7 D n a a 

+ 7 i ^r J i >Z (E i > y x - E i-x^- 


(53) 


p 0 w l a, 2 "" n 0 e n 0 e 


where v, _ was set equal to zero. The first-order equation of motion given by equa- 

A 9 Z 

tion (8) forces v, to vanish. 

In order to compare the amplitudes of the nonlinear fields with those of the linear 
fields, J 2 must first be expressed as a function of Ej and Eg. This is accomplished 
by substituting the expressions for Jj (eqs. (23)) into equation (53), which is the second- 
order Ohm's law. 

Before proceeding, it will be convenient to split the problem into two cases as was 
done for the perpendicular propagation. These two cases arise since there are two 
dispersion relations for the linear wave equation. The wave equation can be expressed as 
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Since is a linear function of Ej, equation (54) represents a set of three homogeneous 

equations. For purely parallel propagating waves Ej x and E l y are not coupled to 

Ej . The dispersion relation for the two possible transverse waves referred to as 
I, z 

right and left circular waves are 


and 



<*> 


(“1 + "c,lK"l ‘ “c.eb 



0 ) 


("1 ’ W c,i )( "l + W c,e ) 


(55) 


(56) 


These expressions are identical to those presented in reference 4. Since both waves are 
independent, the higher-order effects for these two waves will be developed separately. 

The linear currents expressed by equation (23) can be substituted into equation (53) 
since the form of the equations are independent of k (no space derivatives were involved 
in their derivation). Consequently, they can be used with the various modes of propaga- 
tion. When Ei =0, the following expression for the first harmonic current is 
i, z 

obtained: 
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(57a) 


(57b) 


(57c) 


Substituting equations (23), (53), and (57a), (57b), and (57c) into the wave equation 
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(59a) 
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(59b) 
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J 1 ,X J l.x 


+ iE 


l,x 


c,i 


<0 


J l,yri,y 


c,i 


(59c) 


Equation (59c) can be nondimensionalized by multiplying through by Bq/E j and 
expressing the variables in terms of Wj, W , and W g , which results in 


E 2,z E 0 _ K 1 , z W e 
E lx ( W 2-< 


v e i,x yw.i e i, 



+ i + W, 


%y \l 4jtc_ J l, x 

K,i E l,x> 


(60) 


Since (47rc/co .)(J, y /E., ) is invarient to any change in w.,, n, B n , and m. when W, , 

C,1 I) X ^ 2 1 U 1 1 

Wp, and W e remain constant, Eg z Bq/E^ x represents a convenient form in which to 
present the results. 

The relation between kg and k^ is found as was done for the perpendicular wave by 
equating those factors that are spatially and time dependent. This procedure yields equa- 
tion (33). The wave represented by equation (59c) is coupled to the fundamental wave; 
however, the waves represented by equations (59a) and (59b) are not. The terms E 9 

Z, X 

and Eg must vanish since kg is determined from equation (33) and does not satisfy 
equation (59a) and (59b) unless they do. Another possible solution exists for equations 
(59a) and (59b). If the Eg z wave is uncoupled from the Eg wave, a separate wave 
number for Eg ^ can be determined from (59a) and (59b). However, the resulting 
dispersion relation 


k 2_“i 
2 '7 


0) 


( w 2 TW C) i)(¥ w c,e) 


( 61 ) 


is identical to the relation for the fundamental that is assumed present. Consequently, 
no additional information is obtained from (59a) and (59b). The electric field associated 
with the second-order perturbation consequently lies along the magnetic field and is 
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actually a function of the superimposed steady-state magnetic field. This latter charac- 
teristic is unusual since, generally, when a wave is a function of E„, it does not depend 
on the magnetic field. 

In summary, when a parallel propagating linear wave (k parallel to Bq, Ej perpen- 
dicular to Bq) exists in a plasma, the nonlinearity of the equation requires that an addi- 
tional wave be formed whose electric field is perpendicular to that of the linear ordinary 
wave (Eg parallel to Bq). 

As with the perpendicular wave, an amplitude ratio is obtained by dividing the 
second-order perturbation by the square of the amplitude of the fundamental. Hence, 
the amplitude ratio is 


E 


2,z 


i 



( 62 ) 


A nondimensional amplitude ratio l E g z |Bq/|Ej will also be investigated. The 
description of the behavior of the amplitude ratios will be discussed in the section 
RESULTS. 


RESULTS 

In the previous section, general characteristics of the nonlinear wave perturbation 
were noted. The nonlinearity of the equations requires that additional waves are formed. 
For both an ordinary and extraordinary first-order wave, the second-order nonlinear 
disturbance is an extraordinary wave. The magnitudes of higher -order disturbances 
were not examined. 

To determine the relative magnitude of the second-order wave, or the amplitude 
ratio, numerical solutions of the appropriate equations developed in the previous section 
were obtained. The various cases will be discussed individually, but first a few 
general comments are appropriate. For both the ordinary and the right and left circularly 
polarized fundamental, the linear approximation is satisfactory even for large amplitude 
linear components. The disturbances associated with the extraordinary fundamental are 
small above a certain critical frequency; below this frequency, they can be large. 

The quantity of pragmatic interest is not so much the magnitude of the various field 
components but the amplitude of the total electric field. Consequently, the amplitude 
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ratio & of the nonlinear wave is presented in all plots of this report. The nondimen- 

4 

sional amplitude ratio is also presented since it can compactly represent many n and 
Bq combinations. The horizontal bars appearing at the top of the plots represent the 
frequency range where the square of the wave number of the fundamental wave is negative 
and consequently where the assumed periodicity does not exist. 

To obtain the ratio of the amplitudes of the first harmonic to that of the fundamental 
from the figures, 8 is multiplied by |E j | , where |Ej | is the amplitude of the funda- 
mental expressed in statvolts per centimeter. To obtain the amplitude of the first 
harmonic, 8 is multiplied by |Ej | 2 . The regions of the plots of the amplitude ratio 
where nonlinear terms are important will vary with the magnitude of the fundamental. 

Such regions can be defined by the conditions that ^>0.l|E^J^; this condition repre- 
sents a situation where the amplitude of the first harmonic is 1/10 that of the fundamental. 
This condition, however, is where the perturbation technique becomes questionable. 
Consequently, the numerical value of the nonlinear amplitude will be questionable. Even 
though the numerical value is in doubt, these regions represent a condition for which 
nonlinear effects are important. 


Perpendicular Propagation 

Extraordinary wave . - Computer solutions of the amplitude ratio of the extraordinary 

fundamental represented by equation (36) are plotted in figure 1(a) as functions of u>, 

4 11 12 14 1 

for the field intensities of 10 gauss and for densities of 10 , 10 , and 10 centi- 

_3 

meters . (Since it is inconvenient to use nondimensional amplitude ratios in order to 
establish n and B trends, dimensional plots shall be used.) Four singularities are 
present on this plot. The range of influence of all but the first is limited to rather 
restricted frequency regions. At frequencies below the first singularity, the amplitude 

ratio of the first harmonic wave is fairly constant and decreases as the density increases. 

2 3 4 12 

The amplitude ratio for a field strength of 10 , 10 , and 10 gauss and a density of 10 

_ O 

centimeter is presented in figure 1(b). Decreasing the field intensity has shifted the 
first singular point to lower frequencies. Concomitantly, the amplitude ratio for fre- 
quencies less than the first singularity increases significantly as the field decreases. 

The singularity at the lowest frequency occurs when the denominator of the expres- 
sion for the nonlinear field contribution, that is, | E 2 1 (eq. (31)), vanishes. At low 
frequencies, this occurs at the Alfven velocity u> = kc/a^/^ (ref. 1, eq. (3-29)). For 
an atomic hydrogen plasma and for the range of /3's calculated, the first singularity 
occurs at about 0. 2 of the ion cyclotron frequency (fig. 2). At the second singularity 
8 increases because of an increasing kj (k — °°); this corresponds to the "lower 
hybrid" resonance (u> = Wj) of the extraordinary fundamental. The third singularity 
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First-order perturbation frequency, rad/sec 
<b) Charge-particle density, 10^ centimeter^ 
Figure 1. - Concluded. 











Figure 3. - Amplitude ratio g of ordinary wave in atomic hydrogen as 
function of frequency. 1 


occurs when Jg increases without limit. This result occurs when 0)3 « u> c e ; pre- 
sumably there exists a subharmonic of the electron cyclotron resonance. When 
u>i « u> / 2, there exists a nonlinear component with a frequency of u> . This 
component is the one which resonates at . The fourth peak consists of two closely 

Cj 6 

spaced singularities, the first of which occurs at the upper hybrid frequency o> = of 

the extraordinary fundamental. These two poles approached each other as the density 

14 -3 

increased so that at 10 centimeter they could not be numerically resolved. 

These singularities and those found for the other cases suggest that in plasmas 
where particle or wave resonance k = °° exist, such as in experiments of particle 
energy addition or in wave damping in magnetic beaches, nonlinear effects are important. 
The nonlinear effects at these resonances could be realized of course only if dissipative 
phenomena were negligible. This assumption is probably very unrealistic at conditions 
of particle and wave resonance. Consequently, it is expected that the spires appearing 
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in £ would be greatly modified. This conclusion suggests that investigations of dissi - t 
pative wave phenomena should include the nonlinear terms of equations (2) and (3). At 
nonresonant conditions, dissipative effects would have much less influence on wave 
propagation. Consequently, nonlinearity can be evaluated at these conditions with much 
more confidence. The nondimensional amplitude ratio for the extraordinary wave is 

plotted in figure 2. Instead of using W n as a parameter, a proportional quantity defined 

1/2 ” 
as |3 = n ' /Bq is used. 

Ordinary wave . - Computer calculations of the amplitude ratio for the second-order 
perturbation of an ordinary wave (eq. (40)) are plotted in figure 3. As in the previous 
case, this perturbation is also an extraordinary wave. There is, however, an essential 
difference not only in the shape of the curves but more significantly in the magnitude of 
the amplitude ratio. The amplitude of the second-order perturbation originating from an 
ordinary wave is considerably smaller than the amplitude originating from an extra- 
ordinary wave. The amplitude ratio of this nonlinear effect of the ordinary wave is 

3 4 

plotted in figure 3 as a function of frequency for magnetic field intensities of 10 and 10 
gauss and particle densities of 10**, 10'*'/ and 10*^ centimeter - / For the conditions 
plotted, the amplitude ratio is very small compared with that for the case of the extra- 
ordinary wave, except for an extremely narrow frequency band. The half breadth can be 

4 

so narrow that it may not be of any practical interest. As an example, at 10 gauss 
12 -3 

and 10 centimeter , the nonlinear amplitude is greater than the linear component for 
a fundamental of 1 statvolt per centimeter only between the frequencies of 1.05437x10** 
and 1. 05442x10** second - *, a condition that is beyond the accuracy of the model. The 
nondimensional amplitude ratio for the ordinary wave is plotted in figure 4. 


Parallel Propagation 

Calculations of the nonlinear perturbation for the right and left circularly polarized 
parallel propagating wave (eq. (62)) are not presented in figure form since their magni- 
tudes are quite small. Only in the near vicinity of the singularities is equation (62) 
significant; the half breadths are so narrow that they are not of practical interest in the 

present analysis since the breadths could easily be masked by dissipative effects. At 
4 12 3 

10 gauss and 10 centimeter , the nonlinear amplitude is greater than the linear com- 
ponent for a fundamental of 1 statvolt per centimeter only between the frequencies of 
9. 5794521x10"* and 9. 5794566x10"* second - *, a condition which is beyond the accuracy of 
the model. 
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CONCLUSIONS 


If second-order terms are included in the equation of motion and Ohm's law for a 
magnetoplasma and if a small amplitude periodic variation in space and time is assumed, 
waves originate from the nonlinear terms of the equations. These waves are present 
along with the normal linear perturbations in plasmas. The second-order perturbations 
or harmonic waves have the following characteristics: 

1. The first harmonic, if assumed to propagate parallel to k, is an extraordinary 
wave regardless of whether the fundamental is ordinary or extraordinary. In both cases 
it will have the same phase velocity as the fundamental with twice the frequency. 

2. For an extraordinary wave in an atomic hydrogen plasma, the first harmonic of 
the nonlinear effects is characterized by four rather narrow singularities. Another 
region of large nonlinear effects occurs at frequencies lower than the first singularity, 
which occurs at about 0. 2 of the ion cyclotron frequency. The nonlinear effect in this 
region increases when the density increases and/or the magnetic field decreases. Above 
about 0. 2 of the ion cyclotron frequency, the nonlinear effect decreases as the frequency 
increases except in the various frequency regions where singularities occur. 

3. For an ordinary wave, the amplitude of the first harmonic (propagating in the 
direction of the fundamental) is generally small. 

4. For the right and left circularly polarized parallel propagating waves, the non- 
linear effects are negligible except for a very narrow range of frequencies and for a 
large amplitude fundamental. 
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